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ABSTRACT 

Spectral diagrams , illustrating the spectral characteristics of different ujater 
types, mere constructed by means of simple statistical analysis of the various 
reflectance properties of mater areas in southern Scandinavia as registered by 
LANDSAT-1 . There mere indications that mater mhose spectral reproduction is dom- 
inated by chlorophyllous matter (phytoplankton) can be distinguished from mater 
dominated by non-chlorophyllous matter .Differences betmeen lakes, as mell as the 
patchiness of individual lakes , concern! ng secchi disc transparency could be 
visualized after classification and reproduction in black and mhite and in color 
by means of Line Printer, Calcomp Plotter (CRT) and Ink Jet Plotter respectively. 

1 . INTRODUCTION 

Water quality studies mere carried out by densitometer analysis of LANDSAT dia- 
positives mithin different mave length bands at the Remote Sensing Section, 
Department of Physical Geography , Uni versity of Lund in 1974 (Helld^n 1975). 
Thereby it proved possible to trace pollution plumes in the Oresund off Copen- 
hagenhagen and to obtain relative measures of the concentration distribution and 
circulation in the polluted mater under different conditions . The results obtained 
also clearly indicated that there is probably a relation betmeen the grey-tones 
of lakes mithin different spectral bands in LANDSAT imagery and secchi disc 
transparency of the same lakes.lt mas considered valuable for the continued 
mork in this field to test the information bulk in LANDSAT CCT for the appraisal 
of mater qualities by automatic data processing .Parts of Scania and the Oresund 
mere chosen as test areas. The data used for the investigation mere registered by 
LANDSAT-1 on Aug. 8, 1973 (1400-09410). 


2. ERRORS AND THEORY 

According to spectral investigations presented by Fitzgerald (1972) mater pol- 
lutions and inorganic suspensions are characterized by a reflectance maximum in 
0. 6-0.7 microns (MSS 5) and then the reflectance falls rapidly mith increasing 
mave length, mhile the spectral characteristics of phytoplankton are rnore complex. 
Generally speaking chlophyll in mater in the form of plankton absorbs effectiv- 
ely mithin the band 0.4-0. 5 microns mhile reflectance is most intensive in the 
band 0.5-0. 6 microns (MSS 4) mith a secondary minimum in 0.6-0. 7 microns (MSS 5) 
and a secondary maximum in 0.7-0. 8 microns (MSS 6). The pre-requisite for a 
secondary maximum in MSS 6 is that the chlorophyll concentration reaches a lom- 
est critical value in the mater layers closest to the surface since the 
electromagnetic radiation in the near IR is almost entirely absorbed by the 
mater molecules in the uppermost centimeters of the mater body. 
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The radiation that is recorded by the satellite borne sensor system cannot be 
related to the radiation properties of the water exclusively^but the spectral 
properties of the atmosphere must also be taken into consideration . Accordi ng to 
Farrow (1975) about 9^, on an average, of the radiation that reaches the atmosphere 
from the sun is back-scattered into space on a cloudless day. Of course, great 
deviations from this value may occur due to variations in solstice and compo- 
sition of the atmospher e . Provided that 75^ of the radiation emitted from the sun 
reaches a water mass on a cloudless day and that the water re-radiates of 

the radiation, i.e. 7, 5^, via the atmosphere into space, it is evident that the 
atmospheric radiation into space is of great importance since it makes up more 
than half of the radiation registered by the sat elli t e • The atmospheric scat- 
tering, which is partially dependent on wave length (Rayleigh scattering) is most 
intensive in short wave lengths and, then, gradually decreases with increasing 
wave 1 ength , Thus , the atmospheric noise is strongest in MSS 4 where it can some- 
times reach values that are equivalent to the re-radiation of the water, while it 
may be considered neglible in the near IR -bands (MSS 6 - MSS 7)(Farrow 1975, 

NASA 1974). 

In the present investigation no attempt has been made to quantify the effect of 
the atmosphere or solstice on the radiation recorded by LANDSAT-1 since the 
study is based on analyses of one LANDSAT scene only. When studying water quality, 
using LANDSAT data, over long periods and at registrations from different times 
during the year, these factors should be taken into consideration, 

3. DATA AND METHODS 

3.1, DIGITAL PROCESSING OF MULTI SPECTRAL DATA 

Basic research regarding digital processing of pictorial information has been 
carried out at the National Defence Research Institute (FDA) since the 1960s 
and in 1973 this activety was extended to comprise multispectral data. 

The computer processing has been carried out at the Stockholm Computer Center, 
using an IBM 360/75 and also using a PDP-11/40 at FOA. As a necessary com- 
plement to the image processing the FOA computer has special equipment for in- 
put and output of image data, for example three-colour TV, ink jet plotter and 
instruments for converting a transparent film to a digital image or, vice versa. 

FOA have several routines for the conversion of data.These routines are all part 
of the process to convert data of different original types and from different 
storage medium to a form which can be used as a standard routine. 

An important part of the digital processing also consists of manipulations of 
separate spectral bands, e.g. when transforming digital data in order to improve 
the image. Amongst other things it is interesting to investigate the amount of 
picture elements with a certain digital value and to present them in the form of 
histograms or to present pixel by pixel of the digital values registered .Other 
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manipulations luhich may be usad are, for exampl e , wri ti ng only one or a feuj 
selected digital levels u/ith one of the image recorder units available, in- 
cluding line printer , colour coding of the digital levels in a spectral band or, 
in advanced processing to use Fourier transformation technique for the analysis 
of any regular geometric forms in the image . Routi nes for all this manipulations, 
and others, are implemented , or can be included as intergrated parts of the program 
system for the processing of mul tisp ectral digital data at FOA,even if its main 
purpose is simultaneous processing of several ' image layers'. True mul tisp ectral 
processing must also be based on knou/ledge of the statistical characteristics 
of image data-in this case routines are required for the calculations of, for 
exampl e , mean values , variances , co-variances and correlations , generation of 
scattering diagrams etc. .Where generation of images is concerned , colour images 
can , naturally , be created luhere separate spectral bands can be given separate 
colours and also composite images in the form of the difference or the quota 
or more complicated mathematical relations betiueen the pixels uihich correspond 
to the same ground point . Special forms of statistical pre-processing are re- 
quired for automatic classi f ication , i . e . the determination of ujhich "class" the 
pixel belongs to according to rules knou/n in advance and in relation to known 
type objects or in relation to groupings of similar pixels which are determined 
during the course of the calculation. These processing routines are also imple- 
mented in the program system in question. A detailed examination of the available 
processing routines are beyond the scope of this report. Please refer to 
separate reports for the description of this program system and its application 
(Gustafsson and Akersten 1976, Orhaug and Akersten 1976 a, -b, Orhaug et al . 

1976). 

3.2.MFTH0DS OF EVALUATION 

After the selection of subscenes and the definition of extensions by determination 
of corner co-ordinates (line number and pixel), the statistical grey-tone 
distribution for the spectral bands MSS 4 - MSS 7 could be presented in the form 
of histograms and symbol coded plots. The distribution of pixels over digital 
grey- tone levels (DGL) for different types of water was automatically evaluated 

A 

and represented in diagrams .Mean values and standard deviations were calculated 
at the same time. The latter formed the basis of the construction of spectral 
diagrams for the study of spectral characteristics of different water qualities. 

In order to visualize differences in water quality grey-tone maps (line printer 
plots), Calcomp Plots (CRT) and color plots (Ink Jet Plotter) were produced. 

Using the line printer for production of grey-tone plots up to nine different 
grey-tones can be reproduced simultaneously, while up to 15 grey-tone levels can 
be reproduced at the same time in a calcomp plot. The ink jet plotter has the 
capability of producing color images directly on ordinary paper .By forcing ink 
under high pressure through a fine nozzle, small drops of ink are produced. A 
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suitable voltage charges the drops luhich can be modulated by a control electrode. 
With three ink systems having different colours (magenta , yellou; and cyan) color 
images can be produced. The one available at FOA can be fed directly from the 
computer or from a tape drive.lt plots six pixels per millimeter in the colors 
mentioned independently modulateable to 15 intensity levels .Current output format 
allous a 1024 x 1536 pixels image to be plotted on a drum in about 90 seconds. 


4. RESULTS 


4.1. THE COPENHAGEN AREA 

Today most of Copenhagen's waste water is discharged unpurified by two pump 
stations into the Eastern part of Kongedybet at a depth of 6-10 m. The pump 
station situated in the North discharges about 14 mill, m^ waste water per 
annum, while the corresponding quantity for the pump station in the South is 
about 40 mill, m^/year (Dackman et al . 1 971 ) .B esides the sewers mentioned two 
day-water sewers start from Strandvaenget, which, on Aug . 27 , 1 973 , served as sources 
of the waste plume furthest to the North registered by LANDSAT-1 (Helld^n 1975). 
The extension of the plumes and presumed ccncentration distribution are shown in 
Fig. 1. The iso-lines , which are supposed to be a relative measure of the particle 
concentration of the discharge plumes, were drawn with the help of a symbol coded 
line printer plot in MSS 5. The grey-tone distribution in the water areas near 
the coast on all plots reflects not only the varying particle concentration of 
the water but probably also the varying depth of water (Fig. 2-3). 

4. 2. LAKES RINGS30ARNA 

As investigation area and subscene number two the lakes Dstra Rings jbn, V/astra 
Ringsjbn and Satoftasjdn were chosen. 

The differences in water quality between and within the three lakes were visu- 
alized by production of plots of the above mentioned types (Fig. 5-6). The plots 
can also be said to illustrate the extent to which a measurement of the secchi 
disc transparency carried out in situ in any one of the lakes is representative 
of the lake in question. In other words, the applied methodology may be useful for 
obtaining an idea of the patchiness of lakes regarding their secchi disc trans- 
parency. The great variations in reflectance in the different bands, within and 
between the registered lakes, can probably be referred to algae blooming. The 
light strakes , which are most distinct in Lake Ostra Ringsjbn in MSS 6 and MSS 7, 
are interpreted as accumulations of blue-green algae in the most superficial 
water layers . 

4. 3. THE LAKES NW YSTAD 

As investigation area and subscene number three the Lakes Krageholmss jbn , 
Ellestads jbn, Snogeholmss jbn and Sbvdesjbn were chosen. The spectral description 
and varying water quality of the lakes are dealt with under the following 
heading . 
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4. 4. COMPARATIVE SPECTRAL ANALYSES 

Spectral diagrams luere constructed for the respective investigation areas in 
order to carry out comparative studies of the spectral characteristics of 
different water qualities (Fig, 7-8). 

4. 4.1. THE COPENHAGEN AREA 

The spectral functions of the water masses outside Copenhagen are characterized 
by having the highest DGL values in MSS 4 falling with a relatively steep 
gradient to MSS 5. In MSS 5 the gradient becomes more gentle and assumes a per- 
manent value between MSS 5 and MSS 7, The position of the functions in the y- 
direction is supposed to provide a measure of the transparency of the water and 
consequently the digital grey-tone levels , especially in MSS 5, are supposed to 
give a relative measure of the particle concentration in the water. High values 
of DGL correspond to high particle concentrations . The spectral functions con- 
cerning the standard deviations of the waste plumes, on the whole display the 
same characteristics as corresponding functions for DGL mean values. The mutual 
positions of the functions are supposed to reflect the varying homogenity 
(patchiness) of the different water masses regarding transparency and particle 
concentration. 

4. 4. 2. LAKES RINGS30ARNA 

All the functions are characterized by assuming the highest DGL values in MSS 4 
from where they fall to MSS 5 in relatively steep gradients . However , between 
MSS 5 and MSS 6 the gradients decrease and in one case even changes sign, but 
between MSS 6 and MSS 7 they resume values that are equivalent to the course of 
functions between MSS 4 and MSS 5. An embryo of a secondary maximum is thus seen 
in MSS 6 concerning A,B and C, a secondary maximum which is a reality for the 
function D; the latter reflects spectral characteristics of an accumulation of 
superficial plankton (algae) in Lake Ostra Ringsjon.The above mentioned embryo 
or occurrence of a secondary maximum in DGL in MSS 6 is supposed to characterize 
chlorophyllous water. The appearance and mutual position of the functions are 
thereby thought to give information about varying algae concentrations and tran- 
sparencies between the lakes. 

The spectral functions concerning standard deviations are characterized by a 
minimum in MSS 5 and a maximum in MSS 6. The mentioned max. and min. are 
probably characterizing chlorophyllous and, where phytoplankton distribution is 
concerned , heterogenous waters. The mutual positions of the functions reflect 
the varying homogenity , patchiness, of the waters concerning transparency and 
algae concentration .All the functions diverge noticeably from the spectral 
functions of the Copenhagen area. 

4. 4. 3. THE LAKES NW YSTAD 

As in the case of the spectral functions of Lakes Rings jdarna, the characteristics 
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of the lakes (a),(b),(C) and(D)are supposed to be a result of the presence of 
chlorofyll (phytoplankton) in the mater bodies* Homever, the presence of chlor- 
phyll in (b) seems to be almost zero. The presence of chlorophyll gives rise to 
increased reflectance values in MSS 6. Matter not containing chlorophyll in mater 
does not give rise to such an increase (Cf the spectral diagrams of the Copen- 
hagen area). The appearance and mutual positions of the functions are supposed to 
reflect varying algae concentrations and secchi disc transparencies betmeen the 
four lakes. 

The spectral functions regarding standard deviations are, as in the case of Lakes 
Rings jdarna, characterized by a minimum in MSS 5 and a maximum in MSS 6. The min. 
and max. mentioned are supposed to characterize chlorophyllous and, regarding 
phytoplankton distribution, heterogenous maters. 

5. CONCLUSIONS 

The secchi disc transparency of the above-mentioned lakes mas measured in situ 
betmeen Aug. 8 and Aug. 27,1973 by means of a secchi disc, The registrations mere 
carried out by the County Administration of Malmdhus (Nilsson 1974). Since the 
variations in secchi disc transparency in six of the seven lakes only amounted 
to 35 cm and there is little investigation material, it mas not considered approp- 
riate to present any correlations betmeen in situ measurment and DGL.The 
reliability of the secchi disc transparency values can also be questioned as the 
measurements mere taken in only one place in each lake, part of mhich mere appa- 
rently very heterogenous concerning particle concentration and secchi disc trans- 
parency .Consequently it is uncertain to mhat extent the registrations carried out 
in situ mere representative of the respective lakes. 

In spite of the small variations in secchi disc transparency betmeen the lakes in 
question, it is evident from the results obtained that very small differences in 
mater quality can be detected by analysis of LANDSAT digital data. The results 
also suggest possibilities of distinguishing chlorophyllous mater from mater 
mhose secchi disc transparency is above all dependent on suspensions of non- 
chlorophyllous matter. In addition, certain conclusions can be dramn about quality 
conditions at different depths of mater at the same time as a measure of the 
patchiness of different mater masses can be obtained regarding transparency and 
particle concentration. 

In order to determine the statistical relations betmeen different types of 
spectral functions and mater qualities the collection of ground truth data from 
about 200 lakes for comparision mith satellite data is in progress. For further 
investigations laboratory analyses by spectrometer are probably necessary in 
order to determine hom the spectral properties of mater vary mith different 
algae species and mith different concentrations and to mhat extent minerogenic 
suspensions can affect the results. 
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KIG.3, COPENHAGEN WASTE PLUMES. Calcomp plot in MSS 4 - MSS 7, 
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FIG. 5. LAKES R I NGS30AR f\IA , 
generated by Ink Jet 
Plotter in MSS 4 and 
MSS 7 (27/8-73). The 
image is filtered for 
the LANDSAT line scan- 
ner no ise. Green =1 and, 
dark blue=highest 
transparency and loiuest 
partial e(phytoplank ton) 
concentration .Dark red 
=1 oiliest transparency 
and highest particle 
concentration . 
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The left lake is l/astra Ringsjdn,the right lake is Ostra Ringsjdn and the upper- 
most part of Cistra Ringsjon is Lake Satoftasjdn. 
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